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Abstract. We discuss how primordial {e.g. inflationary) non-Gaussianity in the 
cosniological perturbations is left imprinted in the Large-Scale Structure of the 
universe. Our findings show that the information on the primordial non-Gaussianity set 
on super-Hubble scales flows into Post-Newtonian terms, leaving an observable imprint 
in the Large-Scale Structure. Future high-precision measurements of the statistics 
of the dark matter density and peculiar velocity fields will allow to pin down the 
primordial non-Gaussianity, thus representing a tool complementary to studies of the 
Cosmic Microwave Background anisotropies. 
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The most compelling feature of inflation pQ is to provide the seeds for the Large- 
Scale Structure (LSS) of the universe and the anisotropies in the Cosmic Microwave 
Background. Cosmological perturbations are created during inflation from quantum 
fluctuations and "redshifted" to sizes larger than the Hubble radius. The deviation from 
a Gaussian statistics characterizing the cosmological perturbations, i. e. the presence of 
higher-order connected correlation functions, emerges as a key observable to discriminate 
among the various scenarios for the generation of cosmological perturbations produced 
during inflation. Non-Gaussianity (NG) is one of the primary targets of present and 
future Cosmic Microwave Background (CMB) satellite missions j^j. 

The main goals of this Letter are to show how NG - set at primordial epochs on large 
(super- Hubble) scales - propagates to subhorizon scales once cosmological perturbations 
reenter the Hubble radius during the matter- and dark energy-dominated epochs and 
to analyze which signatures primordial NG leaves in the Large-Scale Structure (LSS) 
of the universe through the dark matter density and velocity fields. We start with 
considering the case of flat universe filled by a perfect non-relativistic matter fluid. 
For the purposes of this paper this is a good approximation, although for the scales 
crossing the Hubble radius during the radiation dominated era one should include 
the radiation in the evolution equations, thus using a complete transfer function also 
for those scales. We find it convenient to perform our calculations in the comoving- 
synchronous gauge and then switch to the Poisson gauge [Sj, which allows a more 
direct comparison with the standard Newtonian and Eulerian approach adopted in LSS 
studies. We make use of the formalism developed in Refs. [H Ej which the reader 
is referred to for more details. The synchronous-gauge is defined by the line-element 
ds'^ = a^(r)[— (ir^+7jj(x, r)dx''dx^], where a(r) is the scale factor, r is the conformal time 
and X represent comoving Lagrangian coordinates for the fluid element. The energy- 
momentum tensor is T'^'^ = pu^u", where p is the mass density and, in our comoving 
coordinates, the fluid four- velocity is given by u'^ = (1/a, 0,0,0). The homogeneous 
background (Einstein-de Sitter universe) is described by a scale factor a(r) oc r^. A 
very efficient way to write down Einstein and continuity equations is to introduce the 
peculiar velocity-gradient |I1 "(9* = u\j — = \l^"'l'aji where we have subtracted 
the isotropic Hubble flow. Primes stand for differentiation w.r.t. the conformal time; 
semicolons denote covariant different at ion. From the continuity equation T^'^ .^^ = 0, we 
infer the exact solution for the density contrast 5 = 5p/ p 

5(x,r) = (1 + 5i„(x))[7(x,r)/7in(x)]-V2 - 1 , (1) 

where 7 = det7jj. The subscript "in" denotes the value of quantities evaluated at 
some initial time. Such initial conditions will play a key role once they have been 
properly determined to account for the primordial non-linearities in the cosmological 
perturbations. 

From Eq. it is evident that in the comoving-synchronous gauge the real 
independent degree of freedom is the spatial metric tensor 7jj. The energy constraint 
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reads 

8 24 

^^-^\^'i + -^ + n = —6, (2) 

where TV" ^ is the Ricci tensor associated with the spatial metric with scalar curvature 
TZ = 7^\. The momentum constraint reads -t}^ j^- = -^j; bars stand for a covariant 
differentiation in the three-space with metric 'jij. Finally, the Raychaudhuri equation 

^' + l^ + ^i^:>^ + l^S = 0, (3) 

is obtained from the energy constraint and the trace of the evolution equation 'i?*'^ + 
M'j/t + + - 'i?2)5y4 + TZ'^ - 7^(5*J•/4 = 0. Notice that these equations are 

exact and describe the fully non-linear evolution of cosmological perturbations around 
the Einstein-de Sitter universe (up to caustic formation). In order to show how the 
primordial non-Gaussianities show up in in the matter density contrast, we make a 
perturbative expansion up to second order in the fluctuations of the metric and matter 
variables. The spatial metric tensor can be written as 7,^ = (1 — 2-?/'^^) — ■?/;'^^)) 5ij + 
Xij* + \xfj^i where xlj^ ^"^^ xfj^ are traceless and include scalar, vector and tensor 
(gravitational waves) perturbations. Similarly, we split the density contrast into a linear 
and a second-order part as 5(x, r) = 5'^^)(x, r) + i5'^^)(x, r). At linear order the growing- 
mode solutions for a matter-dominated epoch in the comoving-synchronous gauge are 
given by §\ ^(^)(x,r) = |y^(x) + f^VV(x), and xf^ = (v'.^j ' I'^^j^V) for the 
metric perturbations and S^^^ = ^V^v^ for the linear density contrast. Here V5(x) is the 
so-called peculiar gravitational potential. In writing Xif we have eliminated the residual 
gauge ambiguity of the synchronous gauge as in Ref. P| and we have neglected linear 
vector modes since they are not produced in standard mechanisms for the generation of 
cosmological perturbations (as inflation). We have also neglected linear tensor modes, 
since they play a negligible role in LSS formation. 

Let us now discuss the key issue of the initial conditions which we conveniently fix 
at the time when the cosmological perturbations relevant today for LSS are well outside 
the Hubble radius, i.e. when the (comoving) wavelength A of a given perturbation mode 
is such that A ^ 7Y"^, Ti. = a' /a being the conformal Hubble constant. 

In the standard single-field inflationary model, the first seeds of density fluctuations 
are laid down on super-horizon scales from the fluctuations of a scalar field, the inflaton 
pp. Recently many other scenarios have been proposed as alternative mechanisms to 
generate such primordial seeds. These include, for example, the curvaton and the 
inhomogeneous reheating scenarios where the first density fluctuations are produced 
through the fluctuations of a scalar field other than the inflaton. In order to follow the 
evolution on super-Hubble scales of the density fluctuations coming from the various 
mechanisms, we use the curvature perturbation ( = (^^^ + {l/2)(^'^^ + ■ ■ ■ of uniform 
density hypersurfaces, where C^^'' = —ip^^'' — Ti^y^ and the expression for C*^^) can be 
found in Ref. jSj- The crucial point is that the gauge- invariant curvature perturbation ( 
remains constant on super-Hubble scales after it has been generated duirng a primordial 
epoch and possible isocurvature perturbations are no longer present. Therefore, we may 
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set the initial conditions for LSS at the time when ( becomes constant. In particular, 
,^(2) provides the necessary information about the "primordial" level of NG generated 
either during inflation, as in the standard scenario, or immediately after it, as in the 
curvaton scenario. 

Different scenarios are characterized by different values of C^'^K For example, in 
the standard single-field inflationary model C^^"* = 2 (C^^"*) + ^ (^c ~ 1) HHI, where 
is the spectral index for the scalar perturbations. In general, we may parametrize 
the primordial NG level in terms of the conserved curvature perturbation as in Ref. 
[H], C^^^ = 2ani(C^^^) • In the standard scenario — 1, in the curvaton case 
Oni = (3/4r) — r/2, where r ^ (P(t/p)d is the relative curvaton contribution to the total 
energy density at curvaton decay 2j. In the minimal picture for the inhomogeneous 
reheating scenario, a^i = 1/4. For other scenarios we refer the reader to Ref. [2j. 
The curvature perturbation ( = allows to set the initial conditions for the metric 
and matter perturbations accounting for the primordial contributions. Indeed, in the 
synchronous gauge and on super-Hubble scales the energy density contrast vanishes 
(see, e.g., Ref. ^j) and therefore = 0. From the deflniton of (^^^ it then follows that 
(^(1) ~ — = — 5(y9/3 and at second-order 



where a^i will always signal the presence of primordial NG according to our 
parametrization. One of the best tools to detect or constrain the primordial large- 
scale NG is through the analysis of the CMB anisotropies, for example by studying the 
bispectrum In that case the standard procedure is to introduce the non-linearity 
parameter /ni characterizing NG in the large-scale temperature anisotropies [131 1111121 • 
To give the feeling of the resulting size of /ni when |ani| ^ 1, /ni — — 5ani/3 (see 
Refs. 13 E]; where the sign is opposite for a different convention relating the CMB 
anisotropies to the gravitational potential). 

After having learned how to set initial conditions, our next step is to determine how 
a primordial NG evolves onto subhorizon scales. By perturbing Eq. up to second 
order, we get 



where jij = 5jj+7i| +7i| /2. To compute the metric perturbation ip^'^\ we need only the 
Raychaudhury equation ^ and the energy constraint Perturbing these equations 
up to second oder and inserting the solution for the density contrast (0) in terms of the 




(4) 



+ 1(^(1)^^-^(1) _ (Dj-'^ (1) ) _ i5(i)(^(i)^- _ ^(1)^-) 

' ^ \ / hj Im im jkJ 2 'in j / ' 



(5) 
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The energy constraint also provides the solution for the scalar perturbation of 
the traceless part of the second-order spatial metric xlj i while the vector and 
tensor parts follow by solving the momentum constraint and the remaining evolution 
equation, according to the same procedure outlined in Ref. We find Xif = 

-f{l-ani)T^ipip,ij-f (I - flni) T^ip,iipj + §T^ [(1 - a^Jy^VV + (i - ani)¥''V,fc] % + 

^rV>,fcj + fi(-12^,,,VV + 4(vV)'5^, - f<^''V,fe^5^i) + Ajf, where Ajf, 
describing second-order tensor modes generated by linear scalar perturbations and 
possible time-independent terms arising from the initial conditions, will not be necessary 
for our purposes. Finally, for the density contrast in the synchronous gauge we plug the 
solution (jni) into Eq. Q to obtain 



10 , 
9 



126 

+ (2 - a,i) ^VV] • (7) 

The expression obtained for the metric perturbation ■i/'^^) is characterized by a large- 
scale part (r — > 0) which is dominated by the primordial NG contribution. It is such 
primordial NG which propagates onto smaller scales, once the mode reenters the horizon 
during the matter-dominated epoch. This is evident from the dependence on Oni of the 
term proportional to in Eq. (jU)). 

The last step of our computation consists in expressing the relevant quantities in the 
Poisson gauge i-e. the generalization beyond linear order of the longitudinal gauge, 
by which a more direct comparison with the standard Newtonian approximation adopted 
in LSS observations and N-body simulations is possible. Second-order coordinate 
tranformations may be written as = — ^'f^^^ — {^'^2) ~ ^(1) with = a^'^^ 

and the space shift splits into the scalar and vector parts as = -|- d''^^y The 
density contrast and the velocity perturbation transform correspondingly as in Ref. [Sj. 
The first-order transformation parameters to go from the synchronous to the Poisson 
gauge are given by ^ a^^-* = TLp/3, P^^^ = r'^ip/Q and d^^^"^ = 0. It follows that at linear 
order the metric perturbations are given in terms of the peculiar gravitational potential 
</9 [Sj as V''-^^ = (j)^^^ = (p where (j) = (j)^^^ +(f)^'^^ /2 and 6goo = — a^(l + 2(f)) defines the lapse 
function in the Poisson gauge. Linear tensor perturbations are gauge-invariant while 
the density contrast and the velocity perturbations read [Sj S^^^ = —2(p + r^V^9?/6 and 
v^^^^ = respectively. For simplicity we do not report here the explicit expressions 

we find for the second-order transformation parameters and we give the results for the 
relevant perturbations in the Poisson gauge. We refer the reader to Eqs. (5.12), (5.13) 
and (5.14) of Ref. [S], which give the general formulae for Z^^^-*, df'^ and a'-^-*. We just 
mention here that the initial conditions for the second-order perturbations will show 
up in such a computation. For the density constrast and the peculiar velocity in the 
Poisson gauge we find 
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4 



(2) 

^ " 126 



5(vV)' + 2(^'*V,*i + VvV. 
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/8 \ ^2 36 ■ 
+ 10 (^- - flnij V'VV + y ^ 

(2)i f 6 



+ f (am - f ) - 246 , (8) 
10 /21 



9 V20 + 



+ y y.^'* - 4) , (9) 

where V^G = ^ - |</?'V,i with V^^ = -\[iy^ipf - ip,ikV'% Eqs. ^ and © are the 
main results of this Letter. In particular they clearly show how the primordial NG which 
is initially generated on large scales is transferred to the density contrast and peculiar 
velocity on subhorizon scales. The expression for the density contrast is made of three 
contributions: the standard second-order Newtonian piece (proportional to r^) which is 
insensitive to the non-linearities in the initial conditions; a Post-Newtonian (PN) piece 
(proportional to r^) which carries the most relevant information on primordial NG; the 
super- horizon terms (independent of r). Our findings show in a clear way that the 
information on the primordial NG set on super-Hubble scales flows into the PN terms, 
leaving an observable imprint in the LSS. 

It might also be useful to write the density contrast in Fourier space in terms of the 
linear density contrast, by defining a kernel /C(ki, k2; r), depending on the wavevector 
of the perturbation modes, as 

4(r) = 4' V) + / k2; r)4?C^5(ki + k^ - k) . (10) 



The kernel K, reads 



r;rk, ko- t)-- + 2 (k^J^ 1 (ki • ^2; U2 , ^ 



' ^21,2 2 
^1 2 



- Gfnl(k,,k2;T)jj-^, (11) 



where 



5 63(ki ■ k,) + 172(fc? + kl) 
/ni(ki,k2;r) = - - ani + 



3 ' 42A;2 
{kf + kl){k,.k,) /4 

6 A;2A;2 6 (ki ■ ka^^ 



+ ^ ^\2irr ^^ x(;(k,.k,) + fc^ + fc,^ 



In order to obtain the expression in Eq. (fTTj) we have performed an expansion in 
(^i,2t)^^ ^ 1 up to terms (fcjr)^^ starting from Eq. (jH)). The non-linearity paramater 
/ni defined in this way generalizes the standard definition of Refs. jT3] inferred 
from the Newtonian gravitational potential $ = + /ni (v^^ — (v^^))- Our findings 
indicate that the non-Gaussianity in the density contrast shows a non-trivial shape 
dependence which might help in detecting the PN nonlinearities. In Fig. ^ we plot the 
k-dependence of the non-linearity parameter /ni which can be useful in distinguishing 
the different contributions in the measured bispectrum. The mass-density bispectrum is 
given by {6{k,)6{k2)6{k,)) = (27r)35(3)(ki + ks + k3)[2;C(ki, k^; T)Ps{k,)Ps{k2) + cycl.]. 
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Figure 1. Wavenumber dependence of the non-linearity parameter /„!. The 
continuous and dashed-dotted hnes correspond to the k-dependent part A/ni of /ni 
(/ni for flni = 0), with an angle between ki and k2 of 7r/3 and 27r/3 respectively. The 
horizontal lines show the k-independent part |5ani/3| for a^i — 1, 10,40, 100. 

where Ps{k) is the hnear mass- density power spectrum. By inspecting Eq. (fT^ we 
can roughly estimate the relevance of the primordial NG contribution to the mass- 
density bispectrum (proportional to a^i^l/rk)^) with respect to the Newtonian one. 
For the two to be comparable on a given scale A ~ k~^, reminding the reader that 
/ni ~ — foni (|flni| ^ 1), wc find that the non-linearity parameter /ni must be as large as 
/ni ~ (10^/1 + z) (Mpc/A)^. Thus, we confirm that, on a comoving scale A ~ fewMpc, 
a value of the non-linearity parameter /ni ~ 10^ is required for the primordial NG to 
leave observable effects in the clustering of galaxies at low redshift ^Sl- Going to higher 
redshift would allow to pin down lower values of /ni p^^, provided the galaxy-to-dark 
matter bias evolution is accurately handled PT], or that observables directly related 
to the dark matter density contrast are considered (see e.g. Ref. [H]). As far as the 
velocity divergence is concerned, in the limit of large NG, we also recover the standard 
Newtonian formula relating the NG in the peculiar velocity to the one present in the 
gravitational potential $, V ■ fp = — (2/37^)V^$. Finally, one can easily extend our 
findings for large NG to the case in which dark energy is present. For the density 
constrast, it suffices to replace, in the term proportional to Oni, with AD'^Hq'^Qqj^, 
where is the standard Newtonian linear growth factor of density perturbations ^H] 
and florn is the presenty-day matter density parameter. As for the peculiar velocity, 
one has to replace, in the term proportional to Oni, t with 2D ^TIq'^I-L f {Vtm) ^ where 
fiytj) = d\n D+/d\n a 
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